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Abstract 

We perform fits of unconventional dark energy models to the available data from high- 
redshift supernovae, distant galaxies and baryon oscillations. The models are based either on 
brane cosmologies or on Liouville strings in which a relaxation dark energy is provided by a 
rolling dilaton field (Q-cosmology). An interesting feature of such cosmologies is the possi- 
bility of effective four-dimensional negative-energy dust and/or exotic scaling of dark matter. 
An important constraint that can discriminate among models is the evolution of the Hubble 
parameter as a function of the redshift, H{z). We perform fits using a unifying formula for 
the evolution of H{z), applicable to different models. We find evidence for a negative-energy 
dust at the current era, as well as for exotic-scaling {a~^) contributions to the energy density, 
with 3.3 < 5 < 4.3. The latter could be due to dark matter coupling with the dilaton in 
Q-cosmology models, but it is also compatible with the possibility of dark radiation from a 
brane Universe to the bulk in brane-world scenarios, which could also encompass Q-cosmology 
models. The best-fit model seems to include an a~^-scaling contribution to the energy density 
of the Universe, which is characteristic of the dilaton relaxation in Q-cosmology models, not to 
be confused with the spatial curvature contribution of conventional cosmology. We conclude 
that Q-cosmology fits the data equally well with the ACDM model for a range of parameters 
that are in general expected from theoretical considerations. 
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1 Introduction and summary 



There is a plethora of astrophysical evidence today, from supernovae measurements [1-5] , the 
spectrum of fluctuations in the Cosmic Microwave Background (CMB) [6], baryon oscilla- 
tions [7] and other astrophysical data, indicating that the expansion of the Universe is currently 
accelerating. The energy budget of the Universe seems to be dominated at the present epoch 
by a mysterious dark energy component, but the precise nature of this energy is still unknown. 
Many theoretical models provide possible explanations for the dark energy, ranging from a 
cosmological constant [8] to super-horizon perturbations [9] and time- varying quintessence sce- 
narios [10], in which the dark energy is due to a smoothly varying (scalar) field which dominates 
cosmology in the present era. 

The current astrophysical data arc capable of placing severe constraints on the nature of the 
dark energy, whose equation of state may be determined by means of an appropriate global fit. 
Most of the analyses so far are based on effective four-dimensional Robertson- Walker Universes, 
which satisfy on-shell dynamical equations of motion of the Einstein- Friedman form. Even in 
modern approaches to brane cosmology [11], which arc described by equations that deviate 
during early eras of the Universe from the standard Friedman equation (which is linear in the 
energy density), the underlying dynamics is assumed to be of classical equilibrium (on-shell) 
nature, in the sense that it satisfies a set of equations of motion derived from the appropriate 
minimisation of an effective space-time Lagrangian. 

However, cosmology may not be an entirely classical equilibrium situation. The initial 
Big Bang or other catastrophic cosmic event, which led to the initial rapid expansion of the 
Universe, may have caused a significant departure from classical equilibrium dynamics in the 
early Universe, whose signatures may still be present at later epochs including the present era. 
In this context, there has been proposed a specific model for the cosmological dark energy, 
being associated with a rolling dilaton field that is a remnant of this non-equilibrium phase, 
described by a generic non-critical string theory [12-15]. We call this scenario 'Q-cosmology'. 

Since such a non-equilibrium, non-classical theory is not described by the equations of 
motion derived by extremising an effective space-time Lagrangian, one must use a more general 
formalism to make predictions that can be confronted with the current data. The approach we 
favour is formulated in the context of string/branc theory [16,17], the best candidate theory 
of quantum gravity to date. Our approach is based on non-critical (Liouville) strings [18-20], 
which offer a mathematically consistent way of incorporating time-dependent backgrounds in 
string theory. 

The basic idea behind such non-critical Liouville strings is the following. Usually, in string 
perturbation theory, the target space dynamics is obtained from a stringy cr-model [16] that 
describes the propagation of strings in classical target-space background fields, including the 
space-time metric itself. Consistency of the theory requires conformal invariance on the world 
sheet, in which case the target-space physics is independent of the scale characterising the 
underlying two-dimensional dynamics. These conformal invariance conditions lead to a set 
of target-space equations for the various background fields, which correspond to the Ein- 
stein/matter equations derived from an appropriate low-energy effective action that is invariant 
under general coordinate transformations. Unfortunately, one cannot incorporate in this way 
time-dependent cosmological backgrounds in string theory, since — to low orders in a perturba- 
tive expansion in the Regge slope a' — the conformal invariance condition for the metric field 
would require a Ricci-fiat target-space manifold, whereas a cosmological background necessarily 
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has a non-vanishing Ricci tensor. 

To remedy this defect, and thus be able to describe a time-dependent cosmological back- 
ground in string theory, the authors of Ref. [18] suggested that a non-trivial role should be 
played by a time-dependent dilaton background. This approach leads to strings living in num- 
bers of dimensions different from the customary critical number, and was in fact the first 
physical application of non-critical strings [19]. The approach of Ref. [18] was subsequently 
extended [12-15, 20] to incorporate off-shell quantum effects and non-conformal string back- 
grounds describing other non-equilibrium cosmological situations, including catastrophic cosmic 
events, such as the collision of two brane worlds. 

In a recent work [21], we have presented preliminary constraints on Q-cosmology by means of 
supernova data. An important result of our analysis was the unavoidable presence of negative- 
dust-like contributions to the energy budget of the Universe. Such negative-energy dust terms 
arise naturally in Q-cosmology [21], as a result of the existence of non-equilibrium dark energy 
contributions. However, a similar feature may also characterise certain brane cosmologies as 
a result of effective four- dimensional Kaluza-Klein (KK) graviton mode contributions to the 
energy density of the brane [22]. Such brane models are also characterised by dark radiation 
terms, as a result of non-trivial gravitational bulk dynamics [11]. As discussed in Ref. [21], 
an important distinguishing feature of Q-cosmology is an a~^-scaling contribution to the dark 
energy as a result of dilaton relaxation terms [12]. Such a term is distinct from the spatial 
curvature contribution as we shall review below. The astrophysical data analysis of Ref. [21] 
yields a non-zero value of such dark energy dilaton terms. 

It is the purpose of the current work to update and complete such analyses by taking into 
account the very recent supernova data [4,5], as well as data from (luminous red) galaxies [7, 
23,24] on the evolution of the Hubble parameter as a function of the redshift, which, as we shall 
see, turns out to be an important constraint for model building. In the same manner as in the 
analysis of Ref. [21], we also perform a comparison of the results of the fit for the Q-cosmology 
model with those of the conventional ACDM model and a rival model with super-horizon 
perturbations [9] superposed on an underlying Einstein-Friedman-Robertson- Walker Universe. 
We find that the current data are consistent with both the Q-cosmology and ACDM models. 
On the other hand, for the super-horizon model there appears a 2a incompatibility between 
the best-fit values for the H{z) galactic data and the supernova analysis. As in Ref. [21], the 
current analysis indicates the existence of negative-energy dust and non-trivial dark-energy 
contribution attributed to the dilaton relaxation terms. Our results in this work should be 
considered as complementary to other similar studies of unconventional cosmologies that exist 
in the current literature [25]. 

The structure of the article is as follows. In Sec. [21 the most crucial phenomenological 
aspects of brane and Q-cosmology models are outlined, and a generic parametrization for the 
expansion rate of the Universe, to be used in the experimental fits, is developed, which is 
capable of capturing, in a single formula, the most important features of the various models, 
namely negative-energy dust, exotic scaling components to the energy density and dilaton dark 
energy relaxation terms. The cosmological data analyses are discussed in Sec. [3] for each of the 
three observational sources under study, namely the supernovae, the H{z) measurement from 
distant galaxies observations and the baryon acoustic oscillations. In Sec. HI the astrophysical 
constraints on the cosmological models are combined, leading to the determination of the model 
parameters, and implications on possible unconventional cosmological scenarios are discussed. 
Finally, conclusions and outlook are presented in Sec. [51 
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2 Theoretical models 



2.1 Brane Universes 

We commence our discussion by summarising the basic features of brane cosmology [11]. Ac- 
cording to this picture, matter fields are confined on three-brane worlds, while fields from the 
gravitational multiplet (dilaton, graviton, etc.) are allowed to propagate in the bulk. For defi- 
niteness, in this work we consider five-dimensional models. In such a case, the analogue of the 
effective four-dimensional Friedman equation on the brane reads [11]: 



where pu is the matter density on the brane, and we have identified the brane tension a with 
the Newton's constant SttG/S = cr/18 and the four-dimensional cosmological constant A4 is 
related to the AdS-bulk (negative) cosmological constant A5 as follows: 



For our purposes, we shall assume a fine tuning such that A4 = 0. The term yu/a^ in ([T]) is 
known as dark radiation, and stems from the non-trivial bulk dynamics and energy conservation. 
Quadratic terms in pu do not play a crucial role in late times, and hence can be ignored when 
discussing phenomenology for redshifts below three, which will be of interest to us in this work. 

An interesting feature of brane models has been pointed out in Ref. [22] and concerns 
the contributions to the effective four-dimensional energy density on the brane world by KK 
graviton modes. Such modes appear as massive gravitons on the brane with masses m > 
3H/2, where H represents the expansion rate of a de Sitter brane world we assume here 
for concreteness. The analysis of Ref. [22] has shown that a single very massive KK mode, 
corresponding to a particle with high momentum along the bulk direction, exerts pressure on 
the brane pushing it outwards. As a result, the energy in the bulk decreases leading to a 
decrease in the dark energy term, and therefore an effective negative contribution to the dark 
radiation term. In this way, the KK mode behaves like a negative-energy dust. It should be 
noted at this stage that this analysis is still incomplete and cannot be trivially extended to 
realistic cosmological models, where one should integrate over the entire spectrum of the KK 
modes. Nonetheless, for our phenomenological purposes in this work we assume that KK modes 
behave like negative-energy dust. It should be noted that there is no conflict with positive- 
energy theorems, since the bulk energy density contribution of the KK modes is positive [22]. 

2.2 Q-cosmology 

Next we discuss briefly Q-cosmologies, which is our second class of models to be constrained by 
the data. Details on the dynamical equations describing Q-cosmological models with matter 
and dark energy contributions have been reviewed in Ref. [21], where we refer the interested 
reader. In the analysis that follows, we will only outline the basic predictions for the Hubble 
rate H{z), which will be used in our fits. 

For completeness, however, we stress again that the most important feature of Q-cosmology 
is the off-shellness of the appropriate equations describing the dynamics of the fields in the grav- 
itational string multiplet (graviton and dilatons). The target-space effective action variations 
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with respect to the gravitational field are non zero, and the respective values are determined 
by the requirement that the Liouville mode restores the world-sheet conformal invariance of 
the non-critical string. The identification of the Liouville mode with the target time leads to 
a set of dynamical equations and constraints [12-14], whose consistent solution determines the 
Liouville string Universe, under the standard assumptions of homogeneity and isotropy. 

There are several versions of the Liouville theory, which in general lead to different predic- 
tions. There may be models in which all matter and gravitational fields are non conformal on 
the world sheet, or others in which only the gravitational multiplet deviate from criticality. For 
our purposes in this work and in Ref. [21], we restrict ourselves to the latter category. After the 
identification of the Liouville mode with the target time, which can be explained dynamically 
in such models [26], the pertinent equations for gravitons, for instance near a fixed point in 
theory space (and thus for sufficiently long cosmic times), read: 

X cG, matter 

^ Sg,. = 3^'^ + + 0(^9^) , (3) 

where 5''^'™^**'''' denotes the gravitational and matter effective low-energy action in the target 
space and the dot is a Liouville world-sheet zero-mode derivative, which at the very end is 
identified (dynamically) with the cosmic time. The central charge deficit Q^(t) describes the 
microscopic non-critical model, and is expressed as a functional of graviton g^j^j^ and dilaton $ 
fields of the gravitational string multiplet. To lowest order in the Regge slope a', it is given by: 

Q<Q\t) = ?,g''{(p-{^f + ...), ^^2^-\og^. (4) 

We note at this stage that the positive definiteness of Q^(t) (supercritical string [18]) is a model- 
dependent feature, but it is essential in yielding a time-like Liouville mode, thus allowing its 
interpretation as a target time field. Such supercritical strings may describe, for instance, the 
stringy excitations on the three brane describing the observable Universe in a colliding-brane- 
world scenario [12], with the collision providing the main cause for departure from criticality. 

The off-shell terms of the Eqs. ([3]) imply highly non trivial modifications of the analogue 
of the (effective, four-dimensional) Friedman equation for the Q-cosmology, which now in- 
cludes, apart from the standard matter and dark energy contributions denoted collectively by 
p, also the above-mentioned Liouville-string off-shell corrections, Ap, which are not positive 
definite in general (without affecting, however, the overall positive-energy conditions of this 
non-equilibrium theory) : 

i/^W = ?^p + Ap (5) 

The appropriate parametrisation for H{z) in the Q-cosmology framework at late eras, such 
as the ones pertinent to the supernova and other data (0 < z < 2), where some analytic 
approximations are allowed [21], reads: 

H{z) = Ho {Qsil + z f + Qsil + zY + ^2(1 + zf) , n3 + ns + n2 = l, (6) 

with the densities fl2,3.s corresponding to present-day values {z = 0). It is important to remark 
that, as discussed in detail in Ref. [21], the densities Qi involve mixed contributions from both, 
off-shell dark energy and (dark) matter. In this respect, VL2 denotes the dilaton dark energy 
contribution at late eras, not to be confused with spatial curvature contributions. Similarly, 
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^3 does not denote matter dust, but it includes (negative) dark energy contributions scaling 
with the redshift like dust. This is a result of the coupling of the dilaton as well as higher 
string loop corrections, which, depending on the model, could be significant at late epochs 
of the Universe. In general, there is no conflict with positive energy theorems, in a similar 
spirit to negative-energy dust contributions in brane models, as mentioned previously, which 
can be due to appropriately compactified KK graviton modes [22]. We must note at this 
point that Q-cosmologies can indeed be accommodated in brane-world scenarios. This may 
occur, for instance, in colliding branes, where the initial brane collision constitutes an Early- 
Universe cosmically catastrophic event, responsible for the departure of the string theory on the 
brane from criticality [12]. In such brany Q-cosmologies, KK graviton modes appear naturally, 
contributing negative-energy dust components to the effective four-dimensional energy density 
on the brane. Due to the dilaton dark-energy-relaxation (quintessence) terms, however, in such 
models, one may also obtain in general additional scaling contributions of the type at late 
eras [12,21], which are absent in conventional brane models with constant dilatons (c.f. ([1])). 

Moreover, in Q-cosmologies there are dark matter contributions with exotic scaling, denoted 
by Qs, where the exponent 6 is theoretically close to, but not quite, four; such contributions 
are due to the coupling of dark matter species with the (non-constant) dilaton quintessence 
field. We stress once more that in all the above terms, the off-shell Liouville corrections play 
an important role [14,15]. 

A complete analysis of the non-critical-string and dilaton effects, which turn out to be 
important in the present era after the inclusion of matter, requires a numerical treatment [14,15]. 
The energy budget of the Q-cosmology Universe is at present largely phenomenological, since 
the theoretical analysis, involving string loops, is not sufficiently developed to allow for concrete 
predictions insofar as the values of Qi and 6 are concerned. For instance, in Ref. [15] it was 
assumed, rather ad hoc, that there is no exotic matter present today {z = 0), although this 
is not true in the past. An important constraint is, of course, nucleosynthesis, which in the 
analysis of Ref. [15] has been left essentially undisturbed, occurring at MeV scales. 

In general, in our analysis the three parameters to be determined by the fits are ^3, and S. 
As we shall see, the data seem to point towards the fact that the amount of exotic-scaling dark 
matter today is non zero, in contrast to the assumption of Ref. [15]. This, in turn, implies the 
need for revisiting the phenomenological constraints on supersymmetric particle physics models 
at future colliders, such as the LHC, derived in Ref. [15] within the Q-cosmology framework. 
We shall do so in a future publication. 



2.3 Formulae for data analysis 

In our phenomenological fitting procedure, 6 is left as a free parameter, whose best fit value, 
however, is found, as we shall discuss below, close to four, thus in good agreement with the- 
oretical expectations [13-15]. It is important to notice that such exotic matter scaling has 
important consequences on the relaxation of some of the stringent constraints imposed on su- 
persymmetric models of particle physics from thermal dark matter analyses [15] in standard 
Friedman-Robertson- Walker cosmologies. On the other hand, a value 6 = 4 may correspond to 
dark radiation in brane models [22]. In this sense, our parametrisation ([HD may be considered 
as a unifying formula, which can be used for fitting various Q-cosmologies and brane models: 



H{z) = Ho {n^ii + zf + nsil + zY + n^ii + zf) 



1/2 



(7) 
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with ^3 + fi^ + = 1 and ^3, fi^, 5 as free parameters. 

Before closing this section, we mention that we shall compare our results with two other 
cosmological models: a ACDM cold dark matter model with a cosmological constant [8] , where 
the Hubble parameter is expressed as: 

H{z) = Ho {Quil + zf + + z)3a+-o)) 1/2 ^ 

and the super-horizon model, in which the Universe is assumed to be filled with non-relativistic 
matter only [9] and there is no dark energy of any sort: 

H{z) = [a-^l^ - a^y^^^o) , (9) 

with t) the gravitational potential, \l/£o a free parameter, 1 + z = a~^{t)e^°'^^^~^^^^° and 
a(t)e~"^^*-*''"''*^'' the scale factor of the Robertson- Walker Universe. 

3 Data analysis 

Observational data from three different astrophysical sources — namely supernovae, distant 
galaxies and baryon acoustic oscillations — are analysed in order to set constraints on the afore- 
mentioned cosmological models. The data fits and all analyses involved are performed within 
the ROOT [27] analysis framework and the Minuit [28] minimisation and error computation 
code. 

In what follows, we shall use 'generically' the terminology 'Q-cosmology' whenever we use 
formula ([7]) as a fitting function to refer to both brane models with constant dilatons and 
dilaton-quintessence Q-cosmology models. The reader should bear in mind that the distin- 
guishing features of dilaton-quintessence Q-cosmology models, as compared to brane models 
with constant dilatons, are the a~^-scaling contributions to the energy density of the Universe, 
as well as the exotic dark matter scaling a~^, with 6^4. 

Some important remarks are now in order. In the discussion of ACDM model, we have 
allowed, for the sake of completeness, for non-zero spatial curvature contributions. On the other 
hand, for simplicity, when we discuss Q-cosmology and super-horizon models, we assumed a 
spatially flat Universe. This is reflected in the relevant relations between observables and the 
Hubble parameter to be discussed in the following sections. This should make even clearer the 
distinct role played by the above mentioned dilaton dark energy contributions to H{z), as 
compared to spatial curvature terms of conventional cosmology. 

3.1 Supernovae 

We analyse recent type-la supernovae (SN) data released by the Hubble Space Telescope 
(HST) [4] and the ESSENCE collaboration [5]. Among 16 newly discovered high-redshift 
SNe [4], the so-called 'gold' dataset {Riess07:gold) includes SNe from other sets: 14 SNe discov- 
ered earlier by HST [2], 47 SNe reported by SNLS [3] and 105 SNe detected by ground-based 
discoveries, amounting a total of 182 data points. An additional set of 77 SNe tagged as 'silver' 
due to lower quality of photometric and spectroscopic record is also listed, constituting the 
Riess07:gold+silver dataset. The ESSENCE dataset {WV07), on the other hand, consists out 
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of 60 SNe of 0.015 < z < 1.02 discovered by ESSENCE [5], 57 high-z SNe discovered during 
the first year of SNLS [3] and 45 nearby SNe. Besides these SN sets, a compilation of the 
aforementioned data {WV07+HST) [29], normahsed in order to account for the different hght- 
curve-fitters employed, is eventually analysed. The latter sample, listing 192 SNe, includes all 
high-redshift supernovae (2; > 1) observed so far. 

Supernova data are given in terms of the distance modulus fi = dlogd^ + 25, where the 
luminosity distance di (in megaparsecs) for a spatially flat universe is related to the redshift z 
via the Hubble rate H{z): 

dLiz) = cil + z) (10) 

whereas in the general case of a non-zero spatial curvature contribution the luminosity 
distance is given by [2]: 



sin ^a/|^^a:| /q" j ; foi' a closed Universe 
diiz) = ^^'^/^^ X \ \^Wk\ Jo jjfp)^ for a flat Universe, (11) 

sinh ^a/I^^/cI Jq 5 for open Universe. 



We note that this observable depends on the expansion history of the Universe from z to 
the present epoch, and recall that, although most of the available supernovae have z < 1, 
there is a handful with values of redshift up to ^ ~ 1.8. For illustration purposes, henceforth 
both data and predictions of cosmological models will be expressed as residuals, Afi, from 
the empty-Universe prediction (Milne's model, Qu = 0). The analysis involves minimisation 
of the standard function with respect to the cosmological model parameters, as they are 
introduced via the Hubble expansion rate in Eqs. ([7j), ([8]) and ([9]), for the Q-cosmology, ACDM 
and super-horizon models, respectively. 

The WV07+HST dataset, which amounts to a sample of 192 supernovae in total, is shown 
in Fig. [1] with the respective measurement errors. This sample combines all observed high- 
redshift supernovae, most of them detected by the HST. The predictions of the cosmological 
models under study are also displayed for the best-fit parameter values. 

For the standard ACDM model, two cases are investigated: a spatially fiat Universe, 
where ffTU]) is used to calculate the luminosity distance and the general case of a curved Universe, 
where (fTTj) is employed. The best-fit parameter values, the la errors and the corresponding 

values for the Riess07 (gold and gold+silver) and the WV07+HST datasets are listed in 
Table [H The confidence limits obtained in the (f2M,f^A) plane are shown in Fig. [21 together 
with other astrophysical constraints which will be discussed in the following sections. The fits 
to the three datasets are compatible and in agreement with those presented elsewhere [4,5,29]. 
In both cases investigated — for a spatially fiat and for a curved Universe — the SN data provide 
stringent constraints on the cosmological parameters. 

For the super-horizon model [9], the results of the fits to the three supernova datasets are 
presented in Table [H The favoured parameter values for these datasets are compatible and in 
full agreement with those found in an earlier analysis [21]. The goodness of fit for this model 
is comparable with the one for the standard ACDM model, as shown by comparing the ratio 
of the to the number of degrees of freedom (dof). 
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Figure 1: Residual magnitude of the WV07+HST dataset supernovae [29] versus the redshift. 
Cosmological model predictions are drawn: (i) Milne's Universe (black solid line); (ii) Matter- 
only Universe, VLm = 1 (black dashed- dotted line); (Hi) ACDM model for {Qm,^a) = (0.33,0.85) 
(red long-dashed line); (iv) super-horizon model for ^^^q = —0.90 (blue short-dashed line); and 
(v) Q-cosmology for 6 = 4, = —2.8 and Qg = 0.86 (green thick solid line). 





Flat Universe 


Curved Universe 


Dataset 


^^M x'/dof 


(fiM,^^A) X' xVdof 


Riess07:gold 

Riess07:gold+silver 

WV07+HST 


0.308 ±0.020 160 0.88 
0.295 ±0.017 369 1.43 
0.259 ±0.019 196 1.02 


(0.48,0.95) 156 0.87 
(0.55,1.10) 356 1.39 
(0.33,0.85) 195 1.03 



Table 1: Fits to the ACDM model parameter Qm, assuming a flat Universe, and to parameters 
Qm and VL\, when a curved one is considered. The values favoured by the Riess07:gold, the 
Riess07:gold-I- silver and the WV07+HST datasets are listed, together with the corresponding 
values. 



Dataset 




Riess07:gold 

Riess07:gold±silver 

WV07±HST 


-0.71 ±0.05 167 0.92 
-0.78 ±0.05 383 1.48 
-0.90 ±0.07 200 1.05 



Table 2: Fits to the super-horizon model parameter '^iQ. The values favoured by the 
Riess07:gold, the Riess07:gold-l-silver and the WV07+HST datasets are listed, together with 
the corresponding values. 

In the same manner, Q-cosmology models are compared against the recent supernova data, 
including the off-shell effects, which are parametrised in the Hubble rate ([7]). In general, this 
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Figure 2: Confidence-level contours for the A CDM model for the WV07-hHST sample of 192 su- 
pernovae [29]. The 68.3% C.L. region extracted from the hary on- oscillations observable A [7] 
and the result of the fit to the H{z) measurement [23] from distant galaxies is also superposed. 

parametrisation of the model is described by three parameters, 6, Qs and ^3, which serve as 
free parameters of the fitting function to be determined by the data analysis. The results for 
a simple case where the parameter S is kept fixed and equal to four are listed in Table [31 The 
value (5 = 4 is chosen to match the results of similar previous studies in Q-cosmology [21] for 
the exotic scaling of dark matter contributions, induced by dilatons, and dark radiation in 
brane models without dilatons. The current analysis shows that the Q-cosmology model fits 
equally well the supernova data as the ACDM model. The results for the Riess07:gold and the 
WV07+HST datasets are consistent at the la level. The Riess07:gold+silver dataset yields 
results compatible to the Riess07:gold sample, but slightly different to the WV07+HST. Nev- 
ertheless, this does not constitute a significant discrepancy, since the 'silver' sample comprises 
astrophysical objects whose classification as type-la supernovae is uncertain [4]. The estimated 
parameters are compatible with those derived in previous studies [21]. 



Dataset 


^^3 ^4 ^2 xVdof 


Riess07:gold 

Riess07:gold+silver 

WV07+HST 


-3.4 ±0.8 1.2 ±0.3 3.2 ±0.9 157 0.87 
-4.5 ±0.6 1.61 ±0.25 3.9 ± 0.7 357 1.39 
-2.8 ±0.5 0.86 ±0.22 2.9 ± 0.5 195 1.02 



Table 3: Fits to the Q-cosmology parameters Vt^ and Vt^ for a fixed value 5 = 4. The values 
favoured by the Riess07:gold, the Riess07:gold-hsilver and the WV07-hHST datasets are listed, 
together with the corresponding values. ^2 is determined by the other densities so that 
^4 ± ^3 ± ^2 = 1- 

The analysis was subsequently performed treating 5 as a free parameter. In Fig. [3l the 
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minimum achieved for fixed values of 5 leaving ^3 and Qs unconstrained is plotted versus 
the chosen 6. Despite the fairly low value achieved (x^/dof = 1.029), the goodness of fit does 
not vary strongly with S. For the WV07+HST dataset, the minimum value of is ~ 194.5, 
hence the 68% confidence interval is defined by the < 198 condition, corresponding to a 
quite large range of favoured values for 6. This feature (the weak dependence of the fit results 
on 6) is also evident when the analysis is repeated for the Riess07:gold sample (cf. Fig. [31 
inset), where the best-fit value for 6 shifts from (5^7 (WV07+HST dataset) to 5 ^ 3.1. The 
special case of 5 = 3, yielding a of ~ 181 (Riess07:gold dataset), is of no relevance to the Q- 
cosmology model, corresponding to a generic matter-dominated Universe with no cosmological 
constant, already disfavoured by astrophysical data. To recapitulate, Q-cosmology models fit 
the supernovae data very well, with the analysis yielding a strong correlation between the 
cosmological parameters rather than providing a concrete determination of them. The latter 
may be achieved in conjunction with other astrophysical data, as we shall discuss in Section HJ 
where confidence intervals for this analysis will be presented. 

[Ohc 

CM 

195 
194.9 
194.8 
194.7 
194.6 
194.5 




Figure 3: Minimum and x^/dof for the Q-cosmology model and the WVOl+HST SN dataset 
for various values of the 8 parameter. Inset: the same for the Riess07:gold dataset. 



3.2 H{z) constraints from distant galaxies 

We now attempt to further constrain the cosmological models under study by performing fits 
to the Hubble parameter dependence on redshift, H{z), as determined [23, 24] using mea- 
surements of differential ages of passively evolving galaxies following the method described 
in Ref. [30]. This method was applied on a sample of distant galaxies consisted by the 
Gemini Deep Deep Survey (GDDS) [31] and archival data [32], spanning a redshift range of 
< 2; < 2. The ratio H{z)/Hq of these Hubble-rate measurements to the Hubble constant 
Hq = (73.5 ± 3.2) Mpc^^ kms"^, measured by WMAP3 [6], is shown in Fig. m The point at 
z 0.1 was derived [24] by the Luminous Red Galaxy (LRG) sample reported in the Sloan 
Digital Sky Survey (SDSS) [33]. For comparison, the points derived by the latest supernovae 
data released by the HST [4] are also plotted, even through they were not taken into account 
in the fits discussed in the following. 

Here, the Hubble parameter as a function of the redshift, given in Eqs. (I7j), (IHl) and ([9]) for 
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Figure 4: Reduced Hubble parameter H{z)/Hq as measured by SDSS (square) [24] and distant 
galaxies (triangles) [23] and for the Hq value measured by WMAP3 [6]. The predictions for 
the following cosmological models are drawn: (i) ACDM model for {Qm,^a) = (0.2,0.7) (red 
long-dashed line); (ii) super-horizon model for "if iq = —0.56 (blue short-dashed line); and (Hi) Q- 
cosmology model for 5 = 4, = —1 and fls = 0.25 (green thick solid line). The measurements 
from supernovae data [4] are also shown (circles). 

the three cosmological models, is directly fitted to the aforementioned data points. The model 
predictions for the best-fit parameter values are given in Fig. H] for the three models. 

The ACDM model fits the data very well: if a spatially fiat Universe is assumed, a value 
of Qm = 0.27 ± 0.04 is derived with = 6.47 (x^/dof = 0.809), while in the general case of 
a curved Universe, the best-fit values are flu = 0.2 ± 0.3 and = 0.7 ± 0.6 with = 6.46 
(x^/dof = 0.923). The H{z) curve for the latter parameters is shown in Fig. [5] and the fit 
result with the corresponding Icr-error is superimposed on the SN results in Fig. O It is evident 
that this measurement provides a quite loose constraint on the ACDM model with errors of the 
order of 15% (100%) for a flat (curved) Universe and it merely confirms the supernova analysis 
results. 

On the other hand, the H{z) measurement provides a stringent constraint on the super- 
horizon model parameter: \I'^o = —0.56 ± 0.09 with a fairly good fit of = 6.98 (x^/dof = 
0.872). The respective H{z) curve is drawn in Fig. [51 This \&£o value, however, is more than 
2a higher than the one favoured by the high-redshift supernovae of the WV07-fHST dataset, 
presented in Table [2l An elaborate discussion on this model follows in the next section and in 
Sec. m in view of the baryonic oscillations data. 

For the Q-cosmology model, the same procedure as for the supernovae is followed. For a 
fixed value of 5 = 4, the best-fit parameters are Q3 = 0.2 ± 0.3 and Qs = —1.0 ± 1.0 with 

= 7.57 (x^/dof = 1.08). The fit is not as good as for the other two models, nevertheless it 
is still acceptable with a ratio of x^ over number of degrees of freedom near unity (see also the 
Q-cosmology prediction in Fig. |4]). The errors on the parameters are large, but — as will shall 
see in Sec. S] — they are correlated and the parameters are eventually sufficiently constrained 
especially when combined with other observational sources. 

Likewise, the analysis was repeated treating 5 as a free parameter. The evolution of the 
minimum x^ obtained for various values of 6 (and the corresponding best-fit values for ^3 and 
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fls) is shown in Fig. O The global minimum value of = 7.35 {x^/dof = 1.2) occurs at 5 ~ 3, 
however as for the supernovae, the quality of fit does not change dramatically with 6, the 68% 
confidence level being at = 10.9. This outcome is further discussed in Sec. |4]in conjunction 
with other astrophysical data. 



Q-cosmology - distant galaxies H(z) ] 
8.4 r 




Figure 5: Minimum and x^/dof for the Q-cosmology model and the galactic H{z) determi- 
nation for various values of the 6 parameter. 

The fits to all three models were likewise repeated for the same set of H{z) points with the 
addition of the three points measured by HST. The derived model parameters and the goodness 
of fit were not significantly different than the ones reported above, since the HST points and 
their errors lie along the line connecting the galactic measurements. 

3.3 Baryon acoustic oscillations 

Complementary information on cosmological models is provided by data on Baryon Acoustic 
Oscillations (BAO), which show up in the patterned distribution of galaxies on very large scales 
(> 100 Mpc). The SDSS [7] survey covers the region 0.16 < z < 0.47, i.e. the typical redshift 
is zbao = 0.35. Among other observables, the SDSS collaboration has measured the dilation 
scale, Dy, defined as [7]: 

1/3 

(12) 

where Dm is the co-moving angular diameter distance at redshift z. Equation (fT2l) presupposes 
a single-scale approximation in the redshift-to-distance conversion and an equivalent treatment 
between the line-of-sight and the transverse dilation. These assumptions are valid for the 
fiducial ACDM model for (f^M^^A) — (0.3,0.7) for redshifts around 0.35 [7]. They equally 
apply to the Q-cosmology model under discussion, as the Hubble parameter at 2 ~ 0.35 does 
not deviate substantially from the fiducial ACDM, as shown in Fig. HI The measured value of 
the dilation scale at zbao = 0.35 is ^^(0.35) = 1370 ± 64 Mpc [7]. 

Additionally, the observable ^uh"^, where h is the reduced Hubble constant, has been mea- 
sured by SDSS: ^m^^ = 0.130 ± 0.010. Cosmological models predictions are consequently 
performed by the SDSS collaboration [7] against the observable A, defined as: 

A^D^(.BAO)^^^ (13) 

czbao 
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with a measured value of A = 0.469 ± 0.017. For a constant equation of state of dark energy, 
the quantity A is expressed as [34]: 



A 



^BAO 



^BAO 



nl/3 



\^lk\E{zBAo) 



sin2/3 (v4^/o'^^° ^) > for a closed Universe, 



X < 



2/3 

2BAO dz 
E{z) 

ZBAO dz 
E{z) 



for a flat Universe, 
for an open Universe. 



(14) 



where E{z) = H{z)/Hq is the reduced Hubble parameter. 

The impact of the measurement of observable A on the ACDM model is depicted at la 
level on top of other constraints in Fig. [2l As shown also in Ref. [34], BAO impose a further 
constraint (besides the observations of type-la supernovae) on the ACDM model parameters, 
Qm and If a spatially flat geometry for the Universe is assumed, the matter density is 
determined by observable A to be = 0.273 ± 0.025 [7]. 

In the Q-cosmology case, as mentioned in Sec. [21 the ordinary and dark matter compo- 
nent appears in the Hubble parameter ([7]) mixed with exotic contributions. Therefore, the 
measurement of fluh"^ — and of the observable A consequently — does not constrain the Q- 
cosmology model parameters ^3 and Qs- Nonetheless, it does set limits to the sum of the 
matter constituents of present-day densities fl^ and fig- The predictions of Q-cosmology are 
hence compared against another observable, B, involving the dilation scale Dy reduced so as 
to remove the dependence on the current value of the Hubble constant, Hq: 



B = —Dv{zbao) 
c 



2BAO 



dz 



zbao 



E{z)J E{zbao) 



1/3 



(15) 



For a Hubble constant Hq = (73.5 ± 3.2) Mpc~^kms~^, as measured by WMAP3 [6] — also 
employed in the distant galaxies analysis — , we derive a measured value of B = 0.334 ± 0.021. 
Although, the BAO observations may not provide a determination of the Q-cosmology param- 
eters, nevertheless, as we shall see in Sec. HI they constrain such models when combined with 
other astrophysical data. 

As far as the super- horizon model is concerned, the independent measurement of ^m^^ [7] 
by BAO, implying a matter density Qm — 0.27, practically excludes it, since it predicts a 
matter density equal to unity [9]. Yet the model parameter range allowed by observable B 
(devoid of any dependence on ^m), shown in Fig. [6] (blue curve), marginally agrees with the 
supernovae-favoured value for \E'£o- 



4 Discussion and interpretation 

For the conventional ACDM model, where a redshift-independent equation of state i^de = ~1 
for the dark energy component is assumed, the constraints imposed by supernovae, distant 
galaxies and baryonic oscillations are reviewed in Table [Hand in Fig. [2] as well. The determina- 
tion of the Hubble parameter by high-redshift galaxies loosely constrain the model parameters. 
The other two astrophysical observations, on the other hand, set a combined parameter deter- 
mination around (^m, ^a) — (0.28, 0.75), as shown in Fig. [2l 
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SHCDM model - Baryon oscillations^ 
m 

0.34 



0.32 



0.3 



0.28 
0.26 

-2.5 -2 -1.5 -1 -0.5 

Figure 6: Observable B prediction versus \l/£o for the super-horizon model. The horizonal lines 
show the measured value of B (±la); the left-hand (right-hand) vertical band shows the value 
and ±1(7 error favoured by the supernovae (distant galaxies) data. 





Flat Universe 


Curved Universe 


Data 






Supernovae (WV07+HST) 
Distant galaxies {H{z)) 
Baryon oscillations {A) 


0.259 ±0.019 196 1.02 
0.27 ±0.04 6.47 0.809 
0.273 ±0.025 - 


(0.33,0.85) 195 1.03 
(0.2,0.7) 6.46 0.923 
see Fig. H 



Table 4: ACDM model parameter ^lu, assuming a flat Universe, and parameters one? 
for a curved one resulting from type-la supernovae (WV07+HST), distant galaxies (H{z)) and 
baryon acoustic oscillations (A). 

The results for the super-Hubble model, which is based on cosmological perturbations with- 
out introducing any dark energy component, are summarised in Tableland Fig. [61 The type-la 
supernovae data and the galactic H{z) determination set incompatible constraints at 2a level 
on the model sole parameter, \E'^o- Furthermore, baryon oscillation observations lead to an 
independent measurement of Qm — 0.27, as opposed to the model-predicted unity value for the 
matter density. The "^io values which are consistent with the observable B measurement, on 
the other hand, are barely in agreement with the supernova fit result. In conclusion, the super- 
horizon model appears over-constrained when confronted by different sources of astrophysical 
observations and it is excluded by this analysis. 
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Data 


xVdof 


Supernovae (WV07+HST) 
Distant galaxies {H{z)) 
Baryon oscillations (B) 


-0.90 ±0.07 200 1.05 
-0.56 ±0.09 6.98 0.872 
< -0.8 



Table 5: Super-horizon model parameter "^iq, rendered by type-la supernovae (WV07-I-HST), 
distant galaxies (H{z)) and baryon acoustic oscillations (B). 

We now proceed on discussing the (off-shell) Q-cosmology models. The allowed parameter 
regions are shown in Fig. [7] for various values of the parameter 5 related to the exotic scaling 
of dark matter: 5 = 3.7, 4.1, 4.3. All three astrophysical constraints are presented: the 
WV07±HST supernovae (68.3% and 99.7% C.L.), the distant galaxies (68.3% and 99.7% C.L.) 
and the observable B of BAO (Icr). The (fis,^^^) region for which the parametrisation ([7]) 
for the Hubble parameter is not valid (i.e. H'^{z) < 0) for < 2; < 2 is also indicated. The 
constraints from supernova data largely overlap and lie in parallel with those from baryon 
acoustic oscillations. This stems from the fact that the corresponding observables — luminosity 
distance and dilation scale Dy oc B — are both sensitive to 7^7) • However, the typical 
redshifts z for these measurements are quite different: ^sn = 2 for the supernovae as opposed 
to zbao = 0.35 for BAO, rendering independent bounds. 

The parameter 6 may be restricted by considering the overlap of all three cosmological 
constraints. For values 6 < 3.3, the 68% confidence intervals for BAO and SN do not overlap. 
Furthermore, the la contours for the SN and the galaxies partially coincide only for 6 < 4.3. 
Hence, the allowed range for 6 is fairly wide: 3.3 ^ S ^ 4.3. The approximate values for the 
other two parameters, ^3 and Qs, as defined by the overlapping regions of the la intervals for all 
three astrophysical data are given in Fig.[8]as a function of 6. These values represent the trend of 
the densities Qi with varying 6 rather than a precise determination of model parameters (hence 
the thick curves). To recapitulate, the allowed ranges for the Q-cosmology parameters defined 
in the parametrisation ([7]), as restricted by current type-la supernovae, high-redshift galaxies 
and baryon acoustic oscillations are: 3.3 < 5 < 4.3, —7 < ^^3 < —1.5 and 0.2 < < 4.5 with 
the correlations defined in Fig. [HI 

The negative-energy dust contribution appears essential in order to provide fits to the 
astrophysical data, and, as we have stressed previously, it may not be inconsistent with positive- 
energy theorems of the underlying microscopic cosmological model. In brane-world scenarios, 
for instance, such negative-energy dust might represent effective four- dimensional contributions 
to the energy density of KK-compactification (massive) graviton modes, as measured by an 
observer on the brane, while the corresponding bulk energy-density contributions are perfectly 
positive definite. In Q-cosmology models, such negative energies might also be due to the 
coupling with the non constant dilaton-quintessence field, including higher-string-loop terms, 
again without implying a confiict with positive-energy conditions. In addition to this negative 
dust energy, the data also provide evidence for (positive) exotic-scahng contributions to the 
energy density of the form a~^. The value of 5 = 4 is included in the phenomenologically 
allowed range. Such terms may therefore be interpreted as either dark radiation, in the context 
of brane models [11], or exotic dark-matter contributions in Q-cosmologies [14,15,20], due to 
the coupling with non-trivial quintessence dilatons, or even both types in the case of brany 
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Figure 7: Observational constraints on Q- cosmology for 6 = 3.7, 4.1, 4.3 from top to bottom, 
respectively: (a) 68.3% and 99.7% C.L. contours for the WV07+HST supernova dataset (blue 
solid line); (b) 68.3% and 99.7% C.L. contours for the galactic H{z) (red dashed line); and (c) 
Icr region from the observable B of baryon acoustic oscillations (yellow area). 

Q-cosmology: data-favoured parameter values"! 



-3 
-4 
-5 
-6 
-7 
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Figure 8: Astrophysical data-favoured values of Q-cosmology parameters as defined from the 
overlap of the la regions from supernovae, distant galaxies and baryon oscillations data. 

non-critical-string dilaton-quintessence models [12]. 

Stringent constraints on these brane and off-shell Q-cosmological models may be set by the 
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WMAP measurements [6] of the cosmic microwave background (CMB), especially the position 
of the (acoustic) peaks in the spectrum. Nevertheless, since CMB provides evidence on the last 
scattering surface at a redshift of 2 ~ 1100, more phenomenological input is required in order 
to embark on the pertinent analysis. Reversely, astrophysical data analysis provides a handle 
to exclude regions in parameter space and to guide the phenomenological studies. 

5 Conclusions and outlook 

In this paper, we have performed an analysis on data collected from three distinct astrophysical 
sources, namely type-la supernovae (2; < 2), distant galaxies {z < 2) — via the determination of 
the Hubble expansion rate — and baryon acoustic oscillations observed in luminous red galaxies 
(0.1 ^ z ^ 0.5). We compared them with the respective predictions of the conventional ACDM 
model [8], the super-Hubble model [9] and the Q-cosmology [12-15,21]. The ACDM model, as 
expected, fits the data very well and is mainly constrained by the supernova and BAO data. 
The super-horizon model, on the other hand, is excluded by the BAO-measured matter density 
and yields contradictory (at 2a level) parameter determinations by SN and galactic data. 

Finally, the spatially-flat Q-cosmology predictions, as parametrised in ([7]) for z < 2, fits the 
data very well providing an alternative scenario to account for the dark energy component of 
the Universe. The three model parameters are not univocally determined; their allowed range 
and the correlation between them is rather defined. The parameter 6, associated to the exotic 
scaling of dark matter in the context of Q-cosmology, is restricted in the range 3.3 < 5 < 4.3, 
whilst the present values of (exotic) matter densities, fl^ and fls, vary with S as shown in Fig. [Hi 

The fact that the value of 5 = 4 is included in the allowed range, also points towards 
dark radiation terms in brane models [11]. The data also point towards negative-energy dust 
contributions, which could be due to either dark-energy dilaton terms in Q-cosmologies [12,21], 
or Kaluza-Klein compactification (massive) graviton modes in brane-inspired modes [22]. 

Further detailed studies, such as the theoretical determination of the position of the acous- 
tic peaks in the CMB spectrum, using the underlying formalism of the above models, or the 
measurement of the pertinent (complicated, in general 2;-dependent) equation of state, as, for 
instance, in the analysis of Refs. [23,35], are certainly required in order to provide further 
evidence that might discriminate among the various models. The important point that our 
analysis, however, brought hopefully forward, is the fact that alternative, unconventional Cos- 
mologies are certainly compatible with the current data, and appear on an equal footing to 
the standard ACDM cosmological model. Thus, surprises on the possibility of discovering 
completely new and unconventional physical phenomena in the foreseeable future cannot be 
excluded. 
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